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O ' Abstract 

^ : 

(*^ , The well-known "displace, cut and reflect" method used to generate 

"^ ' disks from given solutions of Einstein field equations is applied to the 

n^. superposition of two extreme Rcissncr-Nordstrom black holes to con- 

struct disks made of charged dust and also non-axisymmetric planar 
{3J[), distributions of charged dust on the z = plane. They are symmet- 

ric with respect to two or one coordinate axes, depending whether 
the black holes have equal or unequal masses, respectively. For these 
^^ ' non-axisymmetric distributions of matter we also study the effective 

potential for geodesic motion of neutral test particles. 
PACS: 04.20.Jb, 04.40.-b, 04.40.Nr 

1 Introduction 

Solutions of Einstein's field equations representing disk-like configurations 
of matter are of great astrophysical interest, since they can be used as mod- 
els of galaxies or accretion disks. Disk solutions can be static or stationary 
and with or without radial pressure. Solutions for static thin disks without 
radial pressure were first studied by Bonnor and Sackfield ^, and Morgan 
and Morgan [2] , and with radial pressure by Morgan and Morgan [2] . Other 
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classes of static thin disk solutions have been obtained [^-[T], while station- 
ary thin disks were studied in |H]-Cni- An exact solution to the problem 
of a rigidly rotating disk of dust in terms of ultraelliptic functions was re- 
ported in jllj . Also thin disks with radial tension jl2| . magnetic fields J13j 
and both electric and magnetic fields mj have been studied. The non-linear 
superposition of a disk and a black hole was first considered by Lemos and 
Letelier |15| . Models of thin disks and thin disks with halos made of perfect 
fluids were considered in ^Hl- The generalization of the "displace, cut and 
reflect" method (Sec. ^ of constructing thick static disks was considered by 
Gonzalez and Letelier J17j . 

An interesting class of solutions of the Einstein-Maxwell field equations 
are the conformastatic space-times with charged dust, in which the charge 
density is equal to the mass density. Therefore, the matter is in equilibrium 
because the mutual gravitational attractions are balanced by the electrical 
repulsions. This kind of matter has been called by some authors "electri- 
cally counterpoised dust" (ECD). Distributions of ECD in equilibrium in 
both classical and relativistic theories were studied by Bonnor JH]. Models 
of oblate and prolate spheroids made of ECD were considered, respectively, 
in ^^ and [^Hl- Giirses J2J and Varela P^I studied static spheres of ECD. 
Disk sources for conformastationary metrics (the stationary version of con- 
formastatic metrics) were considered in jl4| . Although one may intuitively 
expect that astrophysical objects do not have a net charge, there exists 
the possibility that electrons escape from a compact star, leaving behind a 
positively-charged one (see, for example, |23ll2l] ^. 

In this paper we apply the well known "displace, cut and reflect" method 
to the superposition of two extreme Reissner-Nordstrom black holes aligned 
on the z axis to generate static disks of ECD on the z = plane. Next 
we repeat the same procedure to Reissner-Nordstrom black holes aligned on 
the y axis and generate non-axisymmetric distributions of ECD on the plane 
z = 0. We briefly study the effective potential of geodesic motion of neutral 
test particles on these structures. 

The paper is divided as follows. Sec. |21 discusses the Einstein-Maxwell 
equations, the "displace, cut and reflect" method and the particular class of 
conformastat metrics. In Sec. 01 disks of charged dust are constructed using 
the superposition of two aligned Reissner-Nordstrom black holes and in Sec. 
Uwe obtain non-axisymetric distributions of charged dust. We also present 
some analysis of the geodesic motion on these matter distributions. Finally, 
in Sec. El we summarize our results. 



2 Einstein-Maxwell Equations, Disks and Confor- 
mastatic Space- Times 

We consider a static space-time with coordinates {t, x, y, z) and a line ele- 
ment of the form 

ds^ = e'^^^'^'^^dt^ - e^(^'2^'^) (dx^ + dy^ + dz^) . (1) 

The Einstein-Maxwell system of equations is given by 

Gfiu = SttT^u, (2) 



-Lfiu — ^ I ^ fj, ^av + .g^v^ pa^ I ) 



(3) 



F^""-^ = 0, (4) 

where all symbols have their usual meaning. We use geometric units G = 
c = 1. 

The method used to generate the metric of the disk and its material 
content is the well known "displace, cut and reflect" method that was first 
used by Kuzmin ^5j and Toomre [22j to construct Newtonian models of 
disks, and later extended to general relativity (see, for example [SlITU]^. 

The material and electric content of the disk will be described by func- 
tions that are distributions with support on the disk. The method can be 
divided into the following steps: First, in a space wherein we have a com- 
pact source of gravitational field, we choose a surface (in our case, the plane 
z = 0) that divides the space into two pieces: one with no singularities 
or sources and the other with the sources. Then we disregard the part of 
the space with singularities and use the surface to make an inversion of the 
nonsingular part of the space. This results in a space with a singularity 
that is a delta function with support on z = 0. This procedure is math- 
ematically equivalent to making the transformation z -^ \z\ + a, with a 
constant. In the Einstein tensor we have first and second derivatives of z. 
Remembering that dz\z\ = 2'd{z) — 1 and dzz\z\ = 26{z), where 'd{z) and 
6{z) are, respectively, the Heaviside function and the Dirac distribution, the 
Einstein-Maxwell equations give us 

G^, = 87TiT;]r-+Q^u6{z)), (6) 

F^r^^ = 47rr5{z), (7) 



where T^J,™' is the electromagnetic tensor Eq. Q , Q^iv is the energy- momentum 
tensor on the plane z = and J^ is the current density on the plane z = 0. 
For the metric Eq. (^, the non-zero components of Q^y are 

Q\ = ^a^'K. (8) 

Ql = Ql = -^9'%b\ + bl), (9) 

where h^y denotes the jump of the first derivatives of the metric tensor on 
the plane z = 0, 

bpLv = gfjLy,z\z=Q+ ~ 9fiu,z\z=0-^ (10) 

and the other quantities are evaluated at z = 0"*". The electromagnetic 
potential for an electric field is 

A^ = (0,0,0,0). (11) 

Using Eq. Hll|) and Eq. Q, the only non-zero component of the current 
density on the plane z = is 

jt ^ l.f^g^%, (12) 

where a^ denotes the jump of the first derivatives of the electromagnetic 
potential on the plane z = 0, 

dfi = ^fj.,z\z=0+ ~ ^fi,z\z=0--i (13) 

and the other quantities are evaluated at z = 0^. The "physical measure" 
of length in the direction dz for metric (^Q) is ^—gzz^ then the invariant 
distribution is S{z)/ yj—gzz- Thus the "true" surface energy density a and 
pressures or tensions P are: 

'j = v^grzQ\, p = -v^9rzQi = -v^9rzQi. (u) 

Since J'^ = pU^, where U^ = b^ j ^fgu-, the "true" surface charge density p 
is 

P = V-9zzgttJ^- (15) 

Let us now specialize metric Eq. (^ to a conformastat form by consid- 
ering 

e'^i^'V''^ = V-\x,y,z) and 6^^'"'^'''' = V^{x,y,z). (16) 

With this choice, in the absence of matter, the Einstein-Maxwell equations 
are satisfied provided (a) V satisfies Laplace equation and (b) the relation 



between V and the electric potential 0(a;, y, z) is of the form (see, for exam- 
ple, [IZj and the Appendix of |28^ for a detailed deduction) 

(p = ±— + const. (17) 

When static charged dust matter is included, the above conditions imply 
that charge density must be equal to the mass density. An interesting feature 
of conformastat metrics is that they permit the construction of complete 
asymmetrical relativistic configurations of matter in equibrium. With metric 
coefficients given by (J16j) . we find the following expressions for the energy 
density and pressures Eq. ((Tl)) 

P = 0. (19) 

With condition ((T7I), Eq. ^^ gives 



27rV2' 



(20) 



where we have chosen the negative sign in Eq. (|17|) that corresponds to 
positive charges. Thus we have a distribution of charged dust matter where 
the charge density is equal to the mass density on the plane z = 0. We 
use now a particular form of the function V to construct disks and non- 
axisymmetric distributions of ECD. 

3 Aligned extreme Reissner-Nordstrom black holes 
on the z axis 

Electrovacuum solutions representing extreme Reissner-Nordstrom black holes 
in arbitrary positions in static equilibrium have been found by Majumdar 
P^ and Papapetrou 30 . Such solutions can be constructed because, as 
previously stated, the functions V{x, y, z) must satisfy Laplace equation and 
thus can be superposed. Suppose we have two extreme Reissner-Nordstrom 
black holes on z axis, the first with mass mi at z = and the second with 
mass 1712 at z = —Z2 (Fig.^. For convenience we use cylindrical coordinates 
(t, R, z, (/?). The function V representing the superposition is given by 

V = l + ^^^+ , "^^ (21) 

VWT^ y/R^ + {z + Z2)^ 
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Figure 1: Schematic drawing of the configuration used to construct charged 
dust disks. 

If we apply the "displace, cut and reflect" method to the system represented 
by Eq. ((U), we have for the energy density Eq. ((TH|) : 



a 



miaR2 + m2{a + Z2)Tli 



(22) 



where T^i = i?^ + a^ and 7^2 = R'^ + {z2 + a)^- Note that a is always 
non- negative. In Fig. ^ we graph the energy density of the disk Eq. ()22|) 
as a function of R in two situations: in (a) we fix mi = 2, a = 1, ^2 = 1 
and change m2, in (b) we fix mi = 1, m2 = 1, a = 1 and change Z2- In 
situation (a) we note that as m2 is increased, the energy density is more 
uniformly distributed along the radius. In (b) Z2 = corresponds to one 
black hole with mass mi = 2 (curve m2 = in (a)). The energy density 
on i? = assumes a minimum value when Z2 = 1, and then increases as 
Z2 is increased. The curve with Z2 = 1 x 10^ (mass m,2 is very far away) 
corresponds to a disk generated solely by mi. 



4 Aligned extreme Reissner-Nordstrom black holes 
on the y axis 

Now we take two extreme Reissner-Nordstrom black holes, with masses mi 
and m.2, located at {x,y,z) coordinates (0, d/2,0) and (0,— (i/2,0), respec- 
tively (Fig. 131). The function V is written as 



V = l + 



m-i 



+ 



m2 



(23) 



/x2 + (y - ^)2 + 2;2 Jx'^ + {y + iy + z"^ 

Using again the "displace, cut and reflect" method, we obtain an expression 
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Figure 2: Energy density a Eq. I\'22\i as a function of R. (a) mi = 2, a = 1, 
Z2 = 1 are fixed and m2 is clianged from m2 = to rn-2 = 4. (b) mi = 1, 
7712 = 1, a = 1 are fixed and Z2 is changed: Z2 = 0, 1, 3 and 1 x 10^. 




Figure 3: Schematic drawing of the configuration used to construct asym- 
metric charged dust distributions. 



for the energy density Eq. (fTH|l : 



a = 



mi 



R^/VmsR' 



,3/2 



2vr ^R^R^ (^R^ + rni VR^ + m2^/R] 



(24) 



where R 



+ {y 



d/2f + a^ and R2 



X 



2 + (y + d/2f + a2. Eq. (EH) 



represents a non-axisymmetric distribution of matter on the z = plane. In 
general it is symmetric only with respect to the y axis, but when mi = m2 
it is also symmetric with respect to the x axis. The extreme points {xc, yc) 
of the energy density, given by Va = 0, are Xc = and the roots of: 




"^l(yc 



2) 



r5/2 



(25) 



where Ric = {yc — d/2)'^ + a^ and R2C = {yc + d/2)'^ + a^. For ?7ii = m,2 we 
note that 2/c = is always a root of Eq. ((^ . Fig.^a)-(c) is a contour plot 
of Eq. (HH) for mi = 1, m2 = 1, a = 1 and (a) d = 0.5, (b) d = 0.85, (c) 
d = 1. We note that as d is increased, i/c = changes from a maximum to 
a local minimum point, and two other symmetrical extreme points appear. 
This is seen in Fig. E^d), where a is plotted along x = for cases (a), (b) 
and (c). The transition of the extreme point yc = from a local maximum 
to a local minimum is determined by a^yy = evaluated at {xc, yc) = (0, 0). 
In the particular case when ttii = 771-2 = a = k, this becomes 



12k^ 



9d^k^ 



3d^ + ik{k'^ - 2d^y Ak"^ + d^ = 0. 



(26) 
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Figure 4: Level curves of the energy density Eq. H24|). with ttii = 1, m^ = 1, 
a = 1 and (a) d = 0.5, (b) d = 0.85, (c) d = 1. (d) The energy density is 
plotted along x = 0, showing the extreme points yc for cases (a)-(c). 



Eq. (f^ has a real root 



6 



(27) 



Fig. |Sf a)-(c) is another contour plot of the energy density with m-i = 1, 
1712 = 0.5, a = 1 and (a) d = 1, (b) d = 1.3 and (c) d = 1.6. For d ^ 
1.3 there is one maximum extreme point yc- When d ^ 1.3, there appear 
two maximum and one minimum extreme points yc, all in asymmetrical 
positions. This is pictured in Fig. [nijd). 

It is also interesting to study the effective potential of geodesic motion of 
neutral test particles for this kind of matter distribution. The Lagrangean 
associated with the metric (^ is 



2£ = \-^f - Y^{±'^ + f + z^), 



(28) 






(b) 




2 


OQl^ 




1 

y 


(C5 


/70.024 


-1 


v__ 


y 


-2 




■ 



0.04 



0.03 



0.02 



0.01 




Figure 5: Level curves of the energy density Eq. H24|) . with rrii = 1, ?7i2 = 0.5, 
a = 1 and (a) d = 1, (h) d = 1.3, (c) d = 1.6. (d) The energy density is 
plotted along x = 0, showing the extreme points yc for cases (a)-(c). 
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where the dots indicate derivatives with respect to the parameter s. One of 
Lagrange's equations gives 







V^H = E = const. 



rfs \dt) ~ 'dt 
For time-hke geodesic motion on the xy plane, we may write 



y-2p_y2^^2_^y2^ 



•2 , -2 , 1 T-i2 



Thus we can define an efi'ective potential Vcs. as 
For our matter distribution, Eq. H31|) reads 

Veff. = ^ 



1 + 



m-i 



+ 



m2 



(29) 



(30) 



(31) 



(32) 



Calculating V14fr. 
of 



0, one obtains the extreme points Xc = and the roots 



■mi{yc 



+ a' 



3/2 



+ 



m2{yc + 



[yc + 



+ a' 



3/2 



0. 



(33) 



For mi = 7712, yc = is always a root of Eq. (jHSI)- Fig. Eta)-(c) is a contour 
plot of Eq. (inil) for mi = 1, m2 = 1, a = 1 and (a) d = 1, (b) d = 1.4, 
(c) d = 2. As d is increased, yc = changes from a maximum to a local 
minimum point, and two other symmetrical extreme points appear. This is 
seen in Fig. El^d), where V^s, is plotted along x = for cases (a), (b) and 
(c). The transition of the extreme point yc = from a local maximum to a 
local minimum is determined by iyes),yy = evaluated at {xc,yc) = (0,0). 
Setting mi = 7712 = m, we obtain 



(2a^ - d^){4:m + ^/d^^+4a^) = 0, 



(34) 



which has a real root d = a\/2. 

Fig. I7fa)-(c) is another contour plot of Vcff. with mi = 1, 1712 = 0.5, 
a = 1 and (a) d = 2, (h) d = 2.5 and (c) d = 3. For d S 2.5 there is 
one minimum extreme point yc- When d ^ 2.5, there appear two minimum 
and one maximum extreme points yc, all in asymmetrical positions. This is 
pictured in Fig. d^d). 
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Figure 6: Level curves of the effective potential Eq. (|32l) . with m\ = 1, 
7712 = 1, a = 1 and (a) d = 1, (h) d = 1.4, (c) d = 2. (d) The effective 
potential is plotted along x = 0, showing the extreme points yc for cases 
(a)~(c). 
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Figure 7: Level curves of the effective potential Eq. (|5^ . with mi = 1, 
7712 = 0.5, a = 1 and (a) d = 2, (b) d = 2.5, (c) d = 3. (d) The effective 
potential is plotted along x = 0, showing the extreme points yc for cases 
(a)-(c). 
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5 Discussion 

We applied the "displace, cut and reflect" method to a conformastat form of 
the metric to generate disks made of electrically counterpoised dust (ECD) 
in equilibrium. The seed metric used was two extreme Reissner-Nordstrom 
black holes placed along the symmetry axis. Non-axisymmetric matter dis- 
tribution of ECD on the z = plane were obtained from the metric of two 
extreme Reissner-Nordstrom black holes located along the y axis. To our 
knowledge no such configurations of matter satisfying Einstein's field equa- 
tions exist in the literature. We studied the variation of the energy density 
distribution by varying the coordinate distance between the black holes as 
well as their masses. For equal black hole masses the energy density distri- 
bution is symmetric with respect to both the x and y axes, while for unequal 
masses it possesses symmetry only with respect to the y axis. The effective 
potential for geodesic motion of neutral test particles on the plane of such 
a matter distribution was also analysed. 

Of course one could extend the procedure used in this paper to generate 
complete asymmetric planar configurations of ECD by putting N > 2 ex- 
treme Reissner-Nordstrom black holes in arbitrary positions on the xy plane 
and then applying the "displace, cut and reflect" method. 
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